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Abstract 

We calculate the cross sections of the elastic and inelastic breakup 
modes for the inclusive breakup reaction 28 Si( 8 B, 7 Be) at beam energies be- 
tween 10 - 40 MeV/nucleon within a direct fragmentation model formulated 
in the framework of the post form distorted-wave Born-approximation. In 
contrast to the case of the stable isotopes, the inelastic breakup mode is 
found to contribute only up to 30% to the total breakup cross section, which 
is in agreement with the recently measured experimental data. However, 
the high energy tail of the energy spectra of 7 Be fragment is dominated by 
the inelastic breakup mode. The breakup amplitude is found to be domi- 
nated by contributions from distances well beyond the nuclear surface. 
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1 Introduction 



Recently, there has been a lot of interest in the study of the proton drip-line 
nucleus 8 B which is perhaps the most likely candidate for having a proton halo 
structure 0, [|, since its last proton has a binding energy of only 137 keV. 
Several measurements reported lately do seem to provide evidence in favor of this 
possibility. For example, the electric quadrupole moment of 8 B is found to be 
twice as large as the value predicted by the shell model, which can be explained 
with a single particle wave function corresponding to a matter density of root 
mean square (rms) radius of 2.72 fm ||. The observed narrow longitudinal 
momentum (p» ) distribution of the 7 Be fragment emitted in the breakup reaction 
of 1.47 GeV/nucleon 8 B on 12 C target has been interpreted in terms of a greatly 
extended proton distribution in 8 B |j] . The significantly enhanced reaction cross 
sections of 8 B measured at beam energies between 20 - 60 MeV/nucleon are 
shown to be consistent with the large matter radius of 8 B required to explain its 
quadrupole moment 0. 

Nevertheless, the existence of a proton halo in 8 B is still an open issue. Nakada 
and Otsuka have shown that (0 + 2) Two shell model calculations can reproduce 
the measured large quadrupole moment of 8 B without invoking a proton halo 
structure. The interaction cross sections measured by Tanihata et al. f| at 790 
MeV/nucleon are consistent with a normal size of 8 B. The p\\ distribution of the 
7 Be fragment emitted in the breakup reaction of 8 B at 41 MeV/nucleon has been 
found to be dependent on the target mass in Ref. where it is argued that in 
contrast to the situation in the neutron halo nuclei |10| , |i~T|j , the assumption of an 
unusually extended spatial distribution is not necessary to explain the narrow p\\ 
distribution in case of 8 B; the reaction mechanism plays an important role here. 

Breakup reactions in which the halo particle(s) is(are) removed from the pro- 
jectile in the Coulomb and nuclear fields of the target nucleus, have played a 
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significant role in probing the neutron halo structure in some light neutron drip- 
line nuclei (see e.g. for a recent review). The enhanced total Coulomb 



breakup cross sections [|T^, I~4], narrow longitudinal momentum distributions of 
the heavy fragments [H], [15], |lEf , and sharply forward peaked angular distri- 
butions of the valence neutron(s) [I7|, [L8| are some of the pivotal observations 
through which the neutron halo structure has been well manifested. 

Apart from the pn distributions of the 7 Be fragment, some data on the total 
cross section of the breakup reaction 8 B + A — ► 7 Be + X on low mass targets have 
also become available recently j0|. The theoretical studies reported so far have 
used either the Serber type p(J of models JTTJ, |2l| or the diffraction dissociation 
picture ]22| developed by Sitenko and co-workers f23|. Both these approaches are 
essentially semi-classical in nature, hence a more microscopic calculation within 
the quantum mechanical scattering theory is needed to interpret the data prop- 
erly. A proper understanding of the nuclear breakup of 8 B is also important in the 
context of the extraction of the astrophysical S-factor for the radiative capture 
reaction 7 Be(p,7) 8 B from the Coulomb dissociation of 8 B |24], |25| . 

In this paper, we present calculations of the cross sections for the breakup re- 
action 8 B + 28 Si — > 7 Be + X within a direct fragmentation model (DFM), which is 
formulated in the framework of the post form distorted wave Born approximation 



(PFDWBA) |26| , p27| , |28| . As the target nucleus involved in this reaction is very 
light, we shall consider only the nuclear breakup process. However, the Coulomb 
breakup can also be calculated in this theory on the same footing (see e.g. |26| ). 
In the next section we present the details of our formalism. In section 3, our 
results are presented and discussed. The conclusions of our work are described 
in section 4. 
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2 Formalism 



The nuclear breakup cross section consists of two components: the elastic breakup 
(ELB) (also known as "diffraction dissociation") where X corresponds to the 
target nucleus A in its ground state and proton, and inelastic breakup (INELB) 
(also known as " diffraction stripping" ) where X can be any other channel of the 
A + p system. The triple differential cross section for the elastic breakup reaction 
a + A — > c + x + A (e.g for our case, a = 8 B, c = 7 Be, x — p), is defined as 

d 3 CT ELB (a, c) M'aficfJ'x QcQx \^ I nELB |2 /-i \ 

dn c dEjn x (2nnr q a fa 1 p ^ 1 ' 1 ] 

where the transition amplitude fii mi is given by 

xV cx (r cx )u e (r cx )y erne (f cx )x {+) (qa, Ri). (2) 

In Eqs. (1) and (2), i is the orbital angular momentum for relative c + x system 
and yimt are the spherical harmonics. V cx represents the interaction between 
constituents c and x while ut is the wave function for their relative motion in the 
projectile ground state. q a (fi a ), q c (/-tc) and q x (fi x ) are the momenta (reduced 
masses) of the particles a, c and x respectively. x' s denote the scattering wave 
functions which are generated by the appropriate optical potentials in respective 
channels. The system of coordinates used are the same as that given in Ref . [^7j . 



The transition amplitude is a six dimensional integral. By making a zero range 



approximation (ZRA) this integral is reduced to three dimensions [29|, although 
its calculation is still a major problem as it involves a product of three scatter- 
ing waves which converge very slowly. In the ZRA the details of the internal 
structure of the projectile appear in the amplitude only through an overall nor- 
malization constant and the values of i other than zero are necessarily excluded. 
Because of the relative p-state between 7 Be and the proton in the ground state 



of 8 B the ZRA, therefore, is not suitable for this case. However, we introduce a 
constant range approximation (CRA) which reduces the integral in Eq. (2) to 
three dimensions and at the same time allows the non-zero values of i to enter 
in the calculations. We assume that the breakup reaction is strongly peripheral 
and that only those configurations where (1) the proton is in the collinear posi- 
tion between the target nucleus and 7 Be and (2) the relative separation between 
proton and 7 Be is constant (say d fm), contribute to the transition amplitude. 
We can then write 

R 

V(r cx )u e (r cx )y eme (r cx ) = DJ^^- d—)y emt (Ri), (3) 

Hi 

where d is taken to be the distance between origin and the maximum in the 8 B 
ground state wave function. Approximations similar to the CRA have been used 



earlier by Piihlhoffer et al. |3(| and Kubo and Hirata to study the a transfer 
reactions induced by 6 Li and 7 Li projectiles within the DWBA. These authors 
have found a reasonable agreement between the calculated and measured angular 
distributions particularly for momentum matched transitions. Furthermore, as 
long as the reaction is not sensitive to the smaller distances, the results obtained 
with Eq. (3) are in agreement with those of the full finite range calculations. 

Substituting Eq. (3) in Eq. (2), making partial wave expansion for x (<lx, Rx 
and using Eq. (4.13) of [Q, we get 



i x m x LM L 



2L + 1 



x < 4 - mjm e | LM L >< £ x 0£0 | LO > (4) 



where the reduced amplitude (3f^ B is given by 



S B = / dRa<-H<h» aiRi) *'-^' a f l) y L M L {Ri)x {+ \^ Ri). (5) 

In Eq. (5) we have 



and Xe x i s the radial part of the wave function X *{ c Lc-> Ra;)- Eq. (5) is similar to 
the amplitude obtained with the ZRA and can be evaluated by using the method 
described in e.g. pSfl . The expressions for the zero range amplitudes are retrieved 
from Eq. (4) - Eq. (6) by assuming I and d equal to zero. 

The transition amplitude for the inelastic breakup reaction a + A — > c + C, 
where C is some final state of the system A + x, is given by |27j 

y/(2i + 1)PZ? LB = /dr„di2 ( x H *(qc,R/)x a (qL.,B.) 

x V cx (v cx 

where x^lqn, Rx) is the form factor, which is obtained by taking the overlap of 
the wave function for the channel C (which incorporates all the possible reactions 
initiated by the interaction between nuclei A and x), with the wave function 
describing the internal states of the target and projectile nuclei. Using Eq. (3) 
and other steps as described above, an expression similar to Eq. (4) can be 
derived for Pi^ LB with the reduced amplitude given by 

& LB = /^X ( - ) *(qc,aiRi) x£(ga, " 2iZi) yLM,(^)x (+) (qa,Ri), (9) 
J a 2 q x rCi 

where x? x is the radial part of the form factor. Its calculation simplifies greatly 

if we introduce the so called "surface approximation" and write x£ m terms of 

its asymptotic form 

xg{q x ,Ri) = ^S ix>c H^\q x ,R t ), (10) 

where Hj£' = Ge x + tF £x , with F ex and Gi x being the regular and irregular 
Coulomb wave functions. This equation can be rewritten in terms of the elastic 
scattering wave function xt x ( see Eq. (5)) as 

x?Ma>,Ri) = Q Se *' c A xtM,Ri)-FeMRi)), (u) 



where Si Xt i x are the S matrix elements for the elastic channel corresponding to 
the angular momentum l x . The validity of the surface approximation has been 
tested by Kasano and Ichimura who found it to be well fulfilled even for 
the deuteron. We use Eq. (11) for the form factor xf also m the interior region 
in Eq. (9), which is not expected to be a serious approximation as this region 
contributes very little to the whole DWBA integral. In order to calculate the 
inelastic breakup cross section one has to sum over all the channels C 7^ £ x , which 
can be easily done by using the unitarity of the S matrix as all the dependence 
on channel C in the transition amplitude rests in the S matrix Sf Xi c 

e 1 sup r = 1- 1 s, x/x 1 2 (12) 

Thus the inelastic breakup cross section can be written as 

d?a INELB (a, c) _ UaM^xQcQx^/ _REAC lELAS\ I nELB nO |2 /io\ 

dn c dE c <m x - l^^t} ** ' )lPem * ^"'^ 

where /3° is the same as the amplitude defined in Eq. (4) with the wave function 
Xe x replaced by the regular Coulomb function. The partial reaction and elastic 
cross sections crf EAC and &f x LAS are related to the scattering matrix elements 
Si x ,i x by their usual definitions. It may be noted that the quantities required to 
calculate the inelastic breakup are the same as those already calculated in the 
case of elastic breakup. 

In order to study the impact parameter dependence of the breakup cross 
section, we define the "probability of breakup" j*~ u P( a > c ) as 
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/OE(24+l)7t" Ka ' c) (15) 



In Eq. (14), d 3 a(a, c)/dE c dQ c dQ x is the sum of the elastic and inelastic breakup 
cross sections given by Eqs. (1) and (13). 
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3 Results and discussion 



The optical potentials in the entrance and outgoing channels and the constants 
D a and d are required as input in our numerical calculations. Although some 
elastic scattering data for the 8 B, 7 Be + 12 C systems at the beam energy of 
40 MeV/nucleon are available [p4" |, the usual optical model fits to them is un- 
fortunately not reported. Unless otherwise stated, our calculations have been 
performed with the following set of optical potentials, Vr = 123.0 MeV, ry = 
0.75 fm, a v = 0.80 fm, W = 65.0 MeV, r w = 0.78 fm, a w = 0.80 fm, with real and 
imaginary volume Woods-Saxon terms. This potential, which is similar to that 
used in the recent analysis of 9 Li and n Li elastic scattering []35|], has been used 
for both 8 B and 7 Be. The (A 1 ^ 3 + a 1 / 3 ) convention was followed to get the radii 
from the radius parameters. The global Becchetti-Greenlees parameterization 



[ |36|1 was used for the proton-target potential. 

The constants D a (see e.g. JJTJ) and d have been determined with a 8 B ground 
state wave function obtained by assuming it to be a pure OP3/2 proton single 
particle state, with separation energy 0.137 MeV, calculated in a central Woods- 
Saxon potential of geometry, r = 1.54 fm, and a$ = 0.52 fm |J7|]. This gives D a 
= -39.0 MeV fm 3 / 2 and d = 1.8 fm, which has been used in all the calculations 
described in this paper. The radius and diffuseness parameters used by Barker 
[@ ( r o — 1-25 fm, and = 0.65 fm), and Esbensen and Bertsch (r = 1.25 fm, 
and a = 0.52 fm) |T| lead to D a = -40.4 and -41.7 MeV fm 3 / 2 and d = 1.8 
and 1.7 fm respectively. On the other hand, with more elaborate RPA models 
for the 7 Be + p overlap wave function 0, the values of D a and d are found to 
be -38.0 MeV fm 3 / 2 and 1.8 fm respectively. Hence, the constants D a and d do 
not show any marked dependence on the nuclear structure model of 8 B. For the 
sake of comparison with other approaches the simplified potential model as used 
by us seems to be adequate at this stage of the present theory. 
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Fig. 1 shows the calculated elastic (dotted line), inelastic (dashed line) and 
total (solid line) breakup cross sections for the 8 B + 28 Si — > 7 Be + X reaction as 
a function of beam energy together with the data taken from fl9| . We see that 
the measured total breakup cross sections are well reproduced by our calculations 
although the contributions of the elastic and inelastic breakup modes are slightly 
over- and under-predicted respectively. The breakup cross sections decrease with 
beam energy up to 20 MeV/A and after that they are almost constant (although 
the inelastic breakup mode still shows a tendency of decreasing somewhat). The 
nuclear breakup cross sections of 8 B as reported in |E| [19] show a similar type 
of energy dependence in this beam energy regime although their increase below 
20 MeV/nucleon is less pronounced than that seen in Fig. 1. Clearly, more 
measurements are needed to clarify this point. 

A striking feature of the results shown in Fig. 1 is that the contribution of 
the inelastic breakup mode to the total breakup cross section is limited only to 



about 30%, which is in agreement with the experimental data [19|. This is in 
marked contrast to the situation in stable nuclei where this mode of breakup 
makes up about 75 - 80 % of the total (a, c) breakup cross section (see e.g. P8"fl). 
To understand this difference, we show in Fig. 2 the breakup probability (Tg) 
(defined by Eq. (14)) for the reaction 28 Si( 8 B, 7 Be) as a function of the impact 
parameter b (b = (id + 1/2) /q a ). It can be seen that Tg peaks at about 8 fm, 
which is in remarkable agreement with that obtained in Ref. [0 from the semi- 
classical arguments. This is quite large in comparison to the sum of the radii (R s ) 
of 8 B and 28 Si (~ 6 fm). Furthermore, most of the contribution to Tg comes from 
the impact parameters beyond 8 fm, while those from distances < R s are strongly 
suppressed. This clearly shows that the breakup of 8 B takes place far away from 
the nuclear surface which reduces the probability of the inelastic breakup process; 
large impact parameters favor the elastic breakup mode. In contrast, for stable 
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isotopes, the breakup probability peaks around R s (where the inelastic breakup 
mode is maximum) and the drop from the peak value for b > R s is much faster 
than that seen in Fig. 2 |28], |41| . This, explains to some extent the difference in 
the nature of the inelastic breakup cross section of 8 B as compared to that of the 
stable isotopes. The fact that the breakup of 8 B is dominated by contributions 
coming from a large range of impact parameters > R s , is in agreement with 
the observation made earlier in the case of neutron halo nuclei n Li and n Be 
2(| . This could provide an indirect evidence for a larger spatial extension of the 



proton in the ground state of 8 B. It should, however, be stressed that a 0p 3 / 2 
configuration for the p - 7 Be relative motion already leads to a 8 B matter density 
with a larger rms radius J7|, (EJ. Nevertheless, in the present calculations the 8 B 
structure input largely affects only the absolute magnitudes of the cross sections 
(through the constant D a ); the peak position in Fig. 2 is mostly decided by 
the reaction dynamics. It is possible, in principle, to include other components 
(eg. pi/2, iV/2) in the 8 B wave function within this formalism. However, as the 
0^3/2 component carries by far the largest spectroscopic weight (> 90 %), this is 
unlikely to alter our results much. 

As a side remark, we point out that the relative contributions of the elastic 
and inelastic breakup modes are independent of the uncertainties in the values 
of D a and d as the same constants enter in all the cross sections. 

In Fig. 3, we show the contributions of elastic and inelastic breakup modes to 
the energy distribution of the 7 Be fragment at the beam energy of 30 MeV/nucleon. 
One notes that while the elastic breakup mode dominates in the peak region, its 
contribution is very weak towards the high energy end (where the proton energy 
E p — > 0); total cross section is made up mostly of the inelastic breakup mode in 
this region. The threshold behavior of the breakup cross section can be easily 
understood from that of the phase-shift S p of the scattering of the proton from 
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the target. It can be shown that [[41] in the limit E p — > 0, the elastic and inelas- 



tic breakup cross sections are proportional to qp £p+l and R 2e p +l (where R is the 
radius of the nuclear potential) respectively. Therefore, in this limit the elastic 
breakup cross section tends to zero even for the s— wave A + p interaction while 
the inelastic one to a finite value, which explains the observation made in Fig. 3. 
It would be very interesting therefore, to perform measurements for the energy 
spectra of the fragment 7 Be to confirm the E p — > behavior. It may help in fix- 
ing the absolute magnitude of the inelastic breakup cross section in the breakup 
experiments. 

4 Conclusions 

To conclude, we have presented for the first time a fully quantum mechanical 
calculation of the elastic and inelastic modes of the nuclear breakup of 8 B on the 
Si target. We employed a direct fragmentation model which is formulated within 
the framework of the post form distorted wave Born-approximation. This is a 
definite improvement over the semi-classical models of the breakup reactions used 
so far for this purpose. We obtain a good overall description of the experimental 
data measured recently. The inelastic breakup mode is found to contribute only 
up to 30% to the total breakup cross section for the 28 Si( 8 B, 7 Be) reaction, which 
is in contrast to the breakup of the stable isotopes. Most of the contributions to 
the breakup cross section come from the distances far beyond the nuclear surface 
which favors the elastic breakup mode. 

The energy spectra of the 7 Be fragment is dominated by the inelastic breakup 
mode towards the high energy end where the proton energy goes to zero. This 
observation which is beyond the scope of the semi-classical models, is a natural 
outcome of the PFDWBA theory of the breakup reactions and it should be verified 
experimentally. 
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We must stress that the lack of the precise knowledge about the parameters 
of the optical potential for the 8 B, 7 Be + 28 Si system is a potential source of 
uncertainty in our calculations. Therefore, the analysis of the existing data (taken 
at around 40 MeV/nucleon) for this system [53] in terms of the conventional 
optical model would be extremely useful. Moreover, similar studies at other 
beam energies are also clearly needed. The present calculations are not very 
sensitive to the nuclear structure models of 8 B. The direct fragmentation model 
should be improved further so that more detail of the projectile wave function can 
enter in the calculations explicitly. Once such an extended reaction description 
is available, the use of more elaborate nuclear structure models of the projectile 
will be meaningful. It may then become possible to use the breakup data to 
distinguish the wave functions of 8 B obtained from a static potential model from 
those calculated from dynamical approaches || where the deformation of this 
nucleus is taken into account. 

Useful discussions with Hans Geissel is gratefully acknowledged. Thanks are 
also due to Ulrich Mosel for his very kind hospitality to one of us (RS) in the 
University of Giessen. 
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Figure Captions 

Fig. 1. Total cross section for the breakup reaction of 8 B + 28 Si — > 7 Be + X as a 
function of the beam energy. The contributions of the elastic and inelastic 
breakup modes are shown by dotted and dashed lines respectively while 
their sum is depicted by the solid line. The experimental data for the total 
(solid circles), elastic (solid squares) and inelastic (open circles) breakup 
cross sections are taken from |I9| . 

Fig. 2. Breakup probability (T}) (as defined by Eq. (14)) for the reaction studied 
in Fig. 1 at the beam energy of 40 MeV/nucleon as a function of the impact 
parameter. 

Fig. 3. Energy distribution of the 7 Be fragment emitted in the breakup of 8 B on 
28 Si target at the beam energy of 30 MeV/nucleon. The solid, dashed and 
dotted lines have the same meaning as in Fig. 1. 
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Figure 1: Total cross section for the breakup reaction of 8 B + 28 Si — > 7 Be + X 
as a function of the beam energy. The contributions of the elastic and inelastic 
breakup modes are shown by dotted and dashed lines respectively while their sum 
is depicted by the solid line. The experimental data for the total (solid circles), 
elastic (solid squares) and inelastic (open circles) breakup cross sections are taken 
from II . 
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Figure 2: Breakup probability (T}) (as defined by Eq. (14)) for the reaction 
studied in Fig. 1 at the beam energy of 40 MeV/nucleon as a function of the 
impact parameter. 
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Figure 3: Energy distribution of the 7 Be fragment emitted in the breakup of 
8 B on 28 Si target at the beam energy of 30 MeV/nucleon. The solid, dashed and 
dotted lines have the same meaning as in Fig. 1. 
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